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Abstract 

A detailed simulation study is performed for multi-r signatures at the Large Hadron Collider, which can 
be used to probe additional Higgs bosons with lepton-specific Yukawa interactions. Such an extended Higgs 
sector is introduced in some of new physics models at the TeV scale. We here consider the two Higgs doublet 
model with the Type-X Yukawa interaction, where nonstandard Higgs bosons predominantly decay into tau 
leptons. These extra Higgs bosons can be pair produced via s-channel gauge boson mediation at hadron 
colliders; qq ->• Z* -> HA and qq' -> W** -> HH± (AiJ ± ), where H, A and H± are CP -even, odd and 
charged Higgs bosons, respectively. Consequently, multi-r originated signals appear in the final state as a 
promising signature of such a model. We find that the main background can be considerably reduced by 
requiring the high multiplicity of leptons and tau-jets with appropriate kinematical cuts in the final state. 
Thus, assuming the integrated luminosity of a hundred of inverse fb, the excess can be seen in various three- 
and four-lepton channels. With the integrated luminosity of thousands of inverse fb, the determination of 
the mass as well as ratios of leptonic decay branching ratios of these Higgs bosons would also be possible. 
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I. INTRODUCTION 



The SU(3)c x SU{2)i x U{l)y gauge structure of the standard model (SM) for elementary 
particles has been tested precisely |IJ]. A missing piece is the Higgs boson, which is responsible 
for electroweak symmetry breaking and mass generation mechanisms for all SM particles. It is 
expected that the Higgs boson will be discovered at the Large Hadron Collider (LHC) in near 
future. The LHC is also searching for the evidence of new physics beyond the SM [§, g|. The 
LHC has already clarified the absence of light new colored particles; e.g., squarks and gluinos in 
the supersymmetric theories J2|, or forth generation quarks ||. However, a light particle without 
strong interactions has not been ruled out yet by the LHC data because of small production cross 
sections. 

The Higgs sector is totally unknown, since no Higgs boson has been discovered yet In the 

SM, only one scalar iso-doublet field is introduced to spontaneously break the electroweak gauge 
symmetry. However, since there is no reason for the Higgs sector with only one doublet, there 
is a possibility of non-minimal Higgs sectors. There are two important experimental constraints 
on extended Higgs sectors; i.e., the flavor changing neutral current (FCNC) and the electroweak 
rho parameter. In the SM, these constraints are automatically satisfied: FCNC is suppressed by 
the Glashow-Illiopoulos-Maiani mechanism, and the rho parameter is predicted to be unity at the 
tree level due to the custodial SU(2) symmetry. On the other hand, non-minimal Higgs sectors 
suffer from both of them in general. It is known that in the Higgs sector with only doublets, 
the rho parameter is predicted to be unity at the tree level, while Higgs models with higher 
representations such as those with triplet fields predict the rho parameter to be different values 
from unity. Therefore, two Higgs doublet models (THDMs) would be a simplest viable extension 
of the SM. However, in the THDM the most general Yukawa interaction predicts FCNC at the 
tree level, because both the doublet couples to a fermion so that the mass matrix and the Yukawa 
matrix cannot be diagonalized simultaneously. In order to avoid this, a discrete symmetry may 
be introduced under which different properties are assigned to each scalar doublet []?]]. Under this 
symmetry, each fermion couples with only one scalar doublet, and hence there are no FCNC at the 
tree level even in the THDM. 

There are four types of Yukawa interactions depending on the ^-charge assignments; i.e., Type- 
I, II, X and Y. Type-II is the most familiar type of Yukawa interactions in the THDM, which is 
the Higgs sector of the minimal supersymmetric standard model (MSSM). Another interesting 
possibility would be the Type-X THDM, where one Higgs doublet couples with quarks and the 



other does with leptons []8| |1Q|| ■ The Type-X THDM can appear in the Higgs sector of a gauged 



extension of the Type-Ill seesaw model [11|, the model of three- loop seesaw with electroweak 
baryogenesis \12\ and a model for positron cosmic ray anomaly [13|. In the SM-like limit, where 
only one of the CP-even Higgs bosons couples to the gauge bosons, the Yukawa couplings of the 
other Higgs bosons tend to be lepton-specific. Since Yukawa coupling constants are proportional 
to the mass of fermions, these extra Higgs bosons predominantly decay into tau leptons for the 
wide range of the parameter space ||. 

The tau lepton has a relatively short lifetime as compared with the muon. It decays into 



2 



lighter leptons and/or hadrons with neutrinos in the detector. The decay products always produce 
missing energies, which make event reconstructions rather complicated. However, for an energetic 
tau lepton, the missing momentum from its decay tends to be oriented to the same direction of 
the charged track 14]. Therefore, the tau lepton momentum can be approximately reconstructed 
by using the collinear approximation 14]. Furthermore, the decay of the tau lepton is correlated 
by its polarization, which can be used to separate leptonic decays from hadronic decays [15]. 

In this paper, we study multi-r signatures at the LHC in the lepton-specific THDM in the 
SM-like limit. Masses of extra Higgs bosons can be of the order of a hundred GeV under the 
B s Xg'j results. Then the gluon fusion mechanism for such extra Higgs bosons is suppressed 
at hadron colliders. In this case, extra Higgs bosons can be pair produced; qq — > Z* — >• HA and 
qq' — > W^* —7- HH^ (AH^), where H, A and are CP-even, odd and charged Higgs bosons, 
respectively. Produced Higgs bosons mainly decay into tau leptons because the Yukawa coupling 
constant is proportional the fermion mass. These characteristic decay modes can be observed in 
multi-T signatures at the LHC. We perform detailed simulation studies for the pair production of 
the extra Higgs bosons where they subsequently decay into multi-r states. It is found that the 
main background can be considerably reduced by requiring the high multiplicity of leptons and 
tau-jets with appropriate kinematical cuts in the final state. Assuming the integrated luminosity 
of a hundred of inverse fb, the excess can be seen in various three- and four-lepton channels. 

This paper is organized as follows. In Sec. II, we summarize the Type-X THDM and give 
basic constraints on the model. The simulation studies of multi-r signatures at the LHC in the 
lepton-specific THDM are presented in Sec. III. Summary and discussions are given in Sec. IV. 



II. THE MODEL AND CONSTRAINTS 



The Higgs potential of the THDM is defined as [16, 17| 



2 v 1 ' 2 

+ As^iX^) + A^I^X^i) + y [( $ i^2) 2 + (^i) 2 ] , (1) 

where <£j(i = 1,2) are the Higgs doublets with hypercharge Y = 1/2. A softly broken Z<i symmetry 
is imposed in the model to forbid FCNC at the tree level, under which the Higgs doublets are 
transformed as $i — > +3>i and $2 — > —^2 01 • The soft-breaking parameter m| and the coupling 
constant A5 are complex in general. We here take them to be real assuming that CP is conserved 
in the Higgs sector. 

The Higgs doublets can be written in terms of the component fields as 



M Vi+hi 



IZi 



where the vacuum expectation values (VEVs) V\ and V2 satisfy \Jv\ + v\ = v ~ 246 GeV and 
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TABLE I: Variation in charge assignments of the softly broken Z2 symmetry |J . 
tan/3 = V2/V1. The mass eigenstates are obtained by rotating the component fields as 




R(/3) 




R(/3) 




(3) 



where lj^ and z are the Nambu-Goldstone bosons, h, H, A and are respectively two CP-even, 
one CP-odd and charged Higgs bosons, and 



R(0) 



cos 9 — sin ( 
sin 6 cos 9 



(4) 



The eight parameters mf-ml and A1-A5 are replaced by the VEV v, the mixing angles a and 
tan/3, the Higgs boson masses m h ,m H ,m A and m H± , and the soft Z2 breaking parameter M 2 = 
m|/(cos/3sin/3). The coupling constants of the CP-even Higgs bosons with weak gauge bosons 
hVV and HVV(V = W,Z) are proportional to sin(/3 — a) and cos(/3 — a), respectively. When 
sin(/3 — a) = 1, only h couples to the gauge bosons while H decouples. We call this limit as the 
SM-like limit where h behaves as the SM Higgs boson [18, 19 . 

Assuming the discrete symmetry (see TABLE ||) , there can be four types of Yukawa interactions 
in the THDM, i.e., Type-I, II, X and Y g, |; 



/•THDM 

•^yukawa 



Q L Y u $ u u R - Q L Y d <5> d d R - L L Ye$it R + H.c, 



(5) 



where <l?j (/ = u, d or £) is either $1 or $2- In the Type-I THDM, all fermions obtain their masses 
from the VEV of $2- In the Type-II THDM, masses of up- type quarks are generated by $2 while 
those of down-type quarks and charged leptons are acquired by $1. The Higgs sector of the MSSM 
is a special THDM, whose Higgs potential is determined by gauge coupling constants and whose 



Yukawa interaction is of Type-II [16]. In the Type-X Yukawa interaction, all quarks couple to 
$2 while charged leptons couple to $1. Remaining one is referred to the Type-Y. The Yukawa 
interactions are expressed in terms of mass eigenstates of the Higgs bosons as 



C 



THDM 
yukawa 



f=u,d,£ 



(6) 



where Pl(r) are projection operators for left- (right-) handed fermions, and the factors are listed 
in TABLE 0. 
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TABLE II: The mixing factors in each type of Yukawa interactions in Eq. (Q) 



Experimental constraints on masses of Higgs bosons H, A, H^ in THDMs depend on the type 
of the Yukawa interaction. Masses of neutral bosons have been bounded by LEP experiment to be 
m H > 92.8 GeV and m A > 93.4 GeV in the MSSM (Type-II THDM with additional relations) [|. 
At the LHC, neutral Higgs bosons can be produced via gluon fusion gg — > </>° [^] where we define 
(jp = H, A, associated production with heavy quarks pp — > tic/P, bbcjP |2l| ] and the weak boson 
mediated processes pp — > Z* HA H and pp -> W±* -> (jPH* H). For the large tan/3 region, 
stronger mass bounds can be obtained from these production processes at the Tevatron and the 
LHC [23, 24 1 . However, if the Yukawa interactions of H and A are quarkophobic which is realized 
in the wide parameter space in the Type-X THDM, these Higgs bosons are less constrained. The 



search for such Higgs bosons at the LEP experiments is found in Rcf. [25]. On the other hand, 
direct search bounds on charged Higgs boson mass has also been set by LEP experiments as 
m H ± > 79.3 GeV by assuming B{H + — > cs) + B(H + — > t + v) = 1 ^6|. Further stronger bound can 
be obtained in the Type-II(Y) THDM g, |§ from the B s -> X s j results as m H± > 295 GeV (28]. 
The observation of B — > rv decay also constrains the mass of charged Higgs bosons for the large 
tan j3 region p9|]. The LHC can also search for charged Higgs bosons in various production processes 
such as pp 



H + H- |§, gb -»■ tH~ H and pp -> W ± H^ p|. However, these bounds depend 
on the types of the Yukawa interactions. Therefore, the relatively light charged Higgs boson is 
still allowed by experimental data in the Type-I(X) THDM. In the Type-X THDM, the charged 
Higgs boson mass can be constrained by the leptonic decays of tau leptons |33], |34|. Extra Higgs 
boson searches in the MSSM (Type-II THDM) have been well studied so far in the literature. 
Fermiophobic Higgs scenario in the Type-I THDM has also been discussed p5| . 

In this paper, we focus on the Higgs boson search in the Type-X THDM, which is less constrained 
by B decay data. In this model, more than 99% of H and A decay into pairs of tau leptons for 
tan/3 > 3 in the SM-like limit; sin(/3 — a) = 1 The neutral Higgs bosons would be produced 
in pair by qq — > Z* — > HA process at the LHC. These Higgs bosons predominantly decay into a 
four-r state, HA — > (t + t~)(t + t~), which is the characteristic signal of the Type-X THDM. The 
tau leptons further decay into leptons or hadrons with neutrinos. Consequently, there are several 
four-lepton final states such as ££££, £££r^ Wryjh, ^hTyj^ and t^t^t^t^ with missing energies, where 
£ denotes an election e or a muon fi, and is hadronic decay products of the tau lepton. Although 
the branching ratios of the <fi — > decay are not large; namely B((jP — > fi + fi~) ~ (m^/m r ) 2 ~ 

0.35%, the HA — > (fi + h~)(t + t~) decay process would also be a clear signature, since the invariant 
mass of the muon pair has a peak at AL„ ~ m H and Mna — tUa. This signature results in fi + /i~££, 
fi + fi~lTh and fi + fi~ThTh final states. We will give the results for simulations of these signals in the 



5 



next section. 

As for the charged Higgs boson associated production, the processes qq' -> W±* -> (jPH^ would 
be dominant. For tan j3 > 2, more than 99% of charged Higgs bosons decay into a tau lepton with 
a neutrino and 0.35% does into a muon with a neutrino. Therefore, the characteristic signatures 
for this process consist of three tau leptons <fpH — > {t + t~)(t^ 1 v), one tau lepton with two muons 
(j) ^ — > h~)(t ± v) or two tau leptons with one muon (t + t ){ji^v). These signals 

result in the three-lepton final states; HI, llr^, It^t^ and t^t^t/j with the large missing energy. In 
the next section, we study the collider signature of these final states in details. 

III. MULTI-r SIGNATURES 

In this section, we present the results of our simulation study for the multi-r signatures at the 
LHC in order to probe the production of neutral and charged Higgs bosons, which predominantly 
decay into tau leptons and occasionally into muons. 

First, we explain the framework of our simulation and event analysis. Second, we present studies 
for the pair production process of the neutral Higgs bosons. Finally, we present studies for the 
charged Higgs bosons associated production with the neutral Higgs bosons. 



A. Framework for event generation and pre-selection 

The signal events are generated by using MadGraph/MadEvent fj6"|j , where the decay of tau 
leptons is simulated by using TAUOLA [j37T| . The partonic events are passed to PYTHIA |38| for 
parton showering and hadronization. Initial-state-radiation (ISR) and final-state-radiation (FSR) 
effects are included. We choose the collision energy to be 14 TeV, and use the CTEQ6L parton 



distribution functions 3£]. Throughout this paper, we set the masses of extra Higgs bosons to 
m H = 130 GeV, m A = 170 GeV and m H ± = 150 GeV, and take the SM-like limit; sin(/3 - a) = 1. 
The total cross section for pp — > HA is estimated to be 53 fb at the tree level |J. For the charged 
Higgs boson associated production, pp — > c/pH^, the total cross sections are estimated to be 125 fb 
for HH^ production, and 77 fb for AH^ production. These mass splitting among the extra 



Higgs bosons are allowed by electroweak precision data in the SM-like limit |4C, 41]. Background 
events for VV (= ZZ, ZW and WW), ti processes where the weak bosons decay leptonically and 
hadronically, and F+jets (V = W and Z) processes followed by leptonic decays of weak bosons 
are generated by PYTHIA, where the decays of tau leptons are also handled by TAUOLA. The total 
cross sections for these processes are given as 108 pb, 493 pb and 110(30) nb, respectively for VV, 
ti and VK(Z)+jets production processes by PYTHIA. We ignore ET-factor corrections for all the 
signal and background processes, for simplicity. 

First, we perform the pre-selection of the events in various four- and three-lepton channels. In 



6 



order to take into account the detector availability, muons and electrons are required to be isolated 1 
and have p T > 15 GeV and \rj\ < 2.5, where p T is the transverse momentum, the pseudorapidity is 
defined as ij = (1/2) ln[tan(#/2)] from the scattering angle 6 in the laboratory frame. Those muons 
and electrons are counted in the events. Then, we construct primal jets from the final state hadrons 
by anti-£; T algorithm |42| with R = 0.4 using the FastJet package [43]. Among the constructed 



primal jets, we identify the tau-jet candidates by the following criteria; 

a jet with p T > 10 GeV and \r/\ < 2.5 which contains 1 or 3 charged hadrons in a small 
cone (R = 0.15) centered at the jet momentum direction with the transverse energy 
deposit to this small cone more than 95% of the jet. 

The R = 0.4 cone of the primal jet acts as an isolation cone to reduce the mis-tagging probability for 
non-tau jets. We present an estimation of the tau-tagging efficiency and the mis-tagging probability 
in Appendix. |A|. The other jets with p T > 25 GeV and |??| < 5 are regarded as a hard jet, which are 
used to estimate the hard QCD activity of the event. Finally, the missing transverse momentum 
f T is calculated as a negative vector sum of the transverse momentum of visible particles, ft T = 

~ J2vis.PT- 

We assume that electric charges of tau-jets are measurable. Using the charges of tau-jets, 
we require the charge sum of the four leptons in the four-lepton channels to vanish to reduce the 
contributions from the background with mis-identified tau-jets. After the pre-selection of the signal 
and background events, we perform further event analysis for each channel. 



B. Neutral Higgs boson pair production 



Here, we present the results of simulation studies for the neutral Higgs boson pair production 
process. In the Type-X scenario, this process is characterized by the four-lepton signature, where 
the leptons can be charged leptons e, [i and also the tau-jet r^. There are fifteen kinds of the 
four-lepton channels in total. For the convenience of our analysis, we divide the leptonic channels 
into three categories; the channels which contain two or more muons, those which contain two or 
more tau-jets, and those with two or more electrons. However the last channels are difficult to be 
utilized due to the limited statistics and the negligibly small branching ratio of the Higgs bosons 
into electrons. 



1. Four-lepton channels with two or more muons 

As we have mentioned at the end of Sec. ||, the dimuon from the direct decay of the Higgs 
bosons would be a clear signature in the Type-X THDM. Although the decay branching ratio is 



1 The isolation condition for muons is given by Pt^YPt < 0.2 where p5? nc is the sum of the magnitude of the 
transverse momentum of the particles inside the R = 0.4 cone around the muon. The isolation condition for 
electrons is given by Pt/Pt > 0.95, where Pj. is the transverse momentum of the jet which contains the electron 
itself. The jet is constructed from the final state hadrons, electrons, photons and non-isolated muons. For our 
isolation conditions, the finding efficiency of muons is slightly better than that of electrons. 
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FIG. 1: The dimuon invariant mass distribution for the 2/x2tt, channel after the pre-selection (left) and 



after the selection cuts before the m z -window cut in TABLE [I] (right). The solid black line is for the 
signal from HA production, the dashed red line is for the VV production, the dotted green line is for the tt 
production and the dot-dashed blue line is for the y+jets production. In the signal process, m H — 130 GeV 



and m A = 170 GeV are taken. 



only 0.35%, there is no way to ignore this signature. First, we explain the study for the 2fj / 2rf l 
channel in detail, and then we comment and summarize the other channels with two or more muons. 

After the pre-selection, the obtained numbers of events for the 2/i2r^ channel are 91, 441, 393 
and 15675 for the HA, VV, tt and V+jets production processes, respectively, in our simulation 
assuming the integrated luminosity of 100 fb _1 . Notice that at the pre-selection, the charges of 
the muons and the tau-jets are not required to be opposite in each to collect all the 2/i2r^ signals 
through the 4r decay of the Higgs boson pairs. Instead, we require the charge sum of the four 
leptons to vanish. The number of signal events are about two orders of the magnitude smaller than 
background contributions. 

In the left panel in FIG. @, we show distributions of the invariant mass of a muon pair in 
the l^ilTh channel after the pre-selection. The solid black line is for the signal HA production, the 
dashed red line is for the VV production, the dotted green line is for the tt production and the dot- 
dashed blue line is for the V+jets production. Numbers of events for each process are normalized 
so as to correspond to the integrated luminosity of 100 fb -1 . The invariant mass distribution of 
the dimuon for the HA production shows two sharp peaks at ~ 130 GeV and 170 GeV which 
corresponds to a pair of primary muons from the decay of extra neutral Higgs bosons. Due to 
the secondary muon from the decay of tau leptons, the signal events behave a broad distribution 
peaked around 50 GeV as well. The dimuon invariant mass distributions for the VV and V+jets 
processes have a peak at Man — m z through the decay of the Z boson into the dimuon. The 
numbers of events in the distributions for the VV and V+jets processes decrease rapidly at high 
mass regions, but that for the tt process decreases only slowly. 
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FIG. 2: Distributions of kinematical variables for the 2[i2tk channel. Line descriptions follow those in 
FIG. 0. 
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TABLE III: Table for background reductions in the 2fj,2rh channel. Listed are the expected number of events 
for the integrated luminosity of 100 fb -1 at the LHC ^/s = 14 TeV. "O.S." stands for "opposite-sign". 

In FIG. H|, we show the distributions of the transverse momentum pj, for a muon, that pTp for 
a tau-jet, that p^ Th for a tau-jet pair, the missing transverse energy J£t(= \Pt\)i the hadronic 
and leptonic scalar sum of the transverse momentum H^(= X^jet \Pt\) anc ^ H^ p (= Yli Th \Pt\)> 
respectively, in the 2/i2r^ channel after the pre-selection. In all panels in FIG. ^, numbers of the 
signal and background events are normalized to be unity. Line descriptions follow those in FIG. |]. 
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Taking into account the difference of distributions between signal and background events in FIG. ^, 
we employ the following selection cuts; 



p% > 40 GeV, (7a) 

Xt > 30 GeV, (7b) 

< 50 GeV, (7c) 

H^ p > 250 GeV, (7d) 

|M W - m z \ > 10 GeV. (7e) 

Here, we briefly explain the background reduction strategy with these cuts. The transverse mo- 
mentum of tau-jets in signal events is relatively larger than that in the background processes due to 
the heavier mass of Higgs bosons than those of Z and W bosons. Therefore a relatively high p T cut 
for tau-jets could enhance the signal-to-background ratio. We note that the high p T requirement 
of tau-jets is also suitable for the stable tau-tagging efficiency at hadron colliders [44|. Background 



events from the y+jets process contain two mis-identified tau-jets from the ISR jets, with a muon 
pair which comes from the Z/j* — > ji + [i~ decay. Therefore, y+jets background events tend to 
have small 0t, and the cut on^r is expected to reduce the y+jets background significantly. The 
background contribution from the ti events can be reduced by using the cut on because the tt 
events tend to contain many jets due to the b quark fragmentation and ISR/FSR, even though two 
of them are mis-identified as tau-jets. The cut on H^ p can reduce the VV and y+jets backgrounds 
significantly. Furthermore, the events which contain Z — > [i^ [i~ can be reduced by rejecting the 
events with the invariant mass of the muon pair close to m z . 



The results of the signal/background reduction are summarized at each step in TABLE ID. 
We show the expected numbers of events for the integrated luminosity of L = 100 fb _1 for each 
process. In the table, we also include the signal events from the charged Higgs boson associated 
production process pp —> (jpH^. The signal-to-background ratio s/b is evaluated at each step of 
the cuts, where s and b represent the numbers of signal and background events, respectively, taking 
all the HA and (jPH production processes as the signal events. In order to evaluate the signal 



significance, we use the significance estimator S defined as |45] 



S=yj2[(s + b)ki(l + s/b)-s], (8) 

which is also given at each step of the selection cuts. The significance S is proportional to the 
square root of the integrated luminosity. We choose the selection cuts basically to enhance S. 
However, since our background events are estimated based on the leading order cross sections and 
distributions with limited statistics, the s/b ratio should not be small but preferably 0(1) to make 
our results conservative. The largest significance can be obtained after the m^-window cut, where 
the number of the signal events is expected to be about 20 while that of background events is about 
6 giving s/b ~ 3 and S ~ 6 for L = 100 fb _1 . 

In the right panel in FIG. 0, we show the dimuon invariant mass distributions after the selection 
cuts ([Fa|), (|7b[), (7c) and (|7d|). After these cuts, the signal process dominates the total events in the 



2fi2rh channel. Thus the two resonant peaks are more enhanced. The mass resolution is expected 
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to be quite well due to the fine resolution of the muon momentum measurement. On the other 
hand, the expected number of events for L = 100 fb _1 is not sufficient to observe the peaks in the 
distribution. 

Now let us consider the detail of the events for the signal process from HA production. As we 
have mentioned, the 2/i2T/ l final state arises through the HA — > 4r and HA — > 2/i2t decays of the 
Higgs bosons. Assuming the dominant branching ratio of B{(jP — > t + t~) ~ 1, the probability to 
get to the 2fj,2Th final state from the former route is expressed as (2) [B(t — > /i)] 2 [B(r —> t/j)] 2 , 
where the binomial coefficient L) = 6. On the other hand, the probability through the latter route 
is 2B{4>° — > [B{t — > Th)] 2 . Thus, the ratio of the probabilities through the two routes is 

3 [B(t — > fj,)] 2 /B{cfP — > {i + fi~), which is predicted to be about 26 in the Type-X THDM. However, 
in the actual events, the ratio of the expected numbers of the events through the two routes suffers 
from the acceptance cuts on the muons. Especially, the p T cut on the muons is expected to reduce 
the ratio, since the muons from the decay of tau leptons have relatively small transverse momentum. 
From the left panel in FIG. |], the ratio of the number of the excess events at the sharp peaks 
to the number of the events in the continuous distribution is found to be 6.6/84.4 ~ 1./12.8. After 
the selection cuts up to the cut on H^ p in (|7d[), the ratio becomes 1.6/20 ~ 1./12.5 (see the right 
panel in FIG. ||). Thus, we find our selection cuts do not significantly modify the ratio. With this 
fact, the branching ratios of B(<f>° — > can be measured as 

„, , n is Number of the excess events at M UL , peaks r „, , l9 . , 

B(^° -> ^n~) = — — — — ^ v x 3 [B(t ^)] 2 x e, (9) 

JN umber of events m the continuous dist. 

where e is the correction factor which reflects the difference of the kinematical distributions of the 

muons between the 4r and 2^t2r decays of Higgs bosons resulting the difference of the kinematical 

acceptance of the muons. It could be simulated by using the kinematical distributions of the muons 

calculated theoretically and the acceptance cuts for muons. In our simulation, we get e ~ 0.5, thus 

the ratio of the muon acceptance in the 4r decay to that in the 2t2[i decay is V0.5 ~ 0.7. 

Once the peaks in the M,„ distribution are observed with the sufficient number of events, the 

measurement of the branching ratio into the dimuon is possible. If we consider the case where the 

Higgs bosons decay into the other modes as well, such like weak bosons or quarks, the left-hand 

side of Eq. (||) is replaced by the ratio of the branching ratios or the partial decay widths for 

muons and tau leptons, which is exactly the ratio of the square of the Yukawa couplings to muons 

and tau leptons. 

Next, we turn our interest to the tau-jet observables. One of the attractive features of the 2/z2t/ 1 
channel is that the four momenta of tau leptons are reconstructable, if the muon pair comes from 
the direct decay of one of the neutral Higgs bosons. In such a case, the missing momentum in the 
event is expected to come from the hadronic decay of tau leptons, and then the full kinematics 
of tau leptons can be reconstructed by using the collinear approximation for the relation of the 
momenta of the tau leptons and the tau-jets. To enhance the signal events where one of the Higgs 
bosons decays directly into the dimuon, we require the muon pair to have opposite-sign charges. 
Then, we apply the collinear approximation for the tau-jets to determine the four momenta of the 
tau leptons. 
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Here, we briefly explain the collinear approximation to calculate the four momenta of the tau 
leptons. If tau leptons are energetic, the missing momentum from its decay would be along the 
direction of the charged track (either a charged hadron (hadrons) or a charged lepton), p miss ~ 
cp Tj , where p mtss i p T j are the momenta of the neutrino and the charged track, respectively. 
The proportionality constant c can be determined by fixing p miss . Accordingly, the momentum 
of the decaying tau lepton can be approximately reconstructed as p T ~ (1 + c)p T J = z~ 1 p T J , 
where z is the momentum fraction of the charged track from the parent tau lepton. At hadron 
colliders, the transverse components of the missing momentum p T can be measured. Assuming 
that the missing transverse momentum of the event is accounted solely by the missing particles 
in the decays of tau leptons, and applying the collinear approximation for both the tau leptons, 
the missing transverse momentum can be expressed by the momenta of charged tracks, as p T ~ 
p™ ssi + p™ %SS2 ~ ci p?p jl + C2 p T j2 . Unknown parameters c\ and C2 are determined by solving 
simultaneous equations. Using the resulting values of Z\ and Z2, the invariant mass of the tau 
leptons pair is related with that of the tau-jets pair as M% ~ z\ZiM^ T . The fractions Z\ and Z2 
should be between and 1. In our analysis, we set a 10% margin for the upper bound on the cut 
of the momentum fractions, i.e., 0. < z\ 2 < 1.1. This is to take into account the resolution of the 
momentum measurements and also the limitation of the validity of the collinear approximation. 
Actually, by taking the 10% margin the number of the signal events which have the solution is 
enhanced by about 10%, although s/b and S are almost unchanged. The lower limit of z is not 
significant because of the high p T cut on the tau-jets. 

In the left panel in FIG. ||, we show the M ThTh distributions, after requiring opposite-sign for 



the muon pair in TABLE [Hj. In the right panel in FIG. |3j, the M TT distributions of the pair of the 
reconstructed tau leptons from the two tau-jets are shown, while the events without the solution 
of the equations in the < < 1-1 ranges are rejected. Similarly to the M MM distribution, two 
resonant peaks can be observed with relatively wide widths at the masses of the neutral Higgs 
bosons. Only a few backgrounds from VV and tt production are expected. The VV background 
can be further reduced by the cut on the m^-window for the reconstructed M TT . 

In the upper panels of FIG. ||], the double invariant mass distribution of M ThTh vs. (left) 
and M TT vs. Muu (right) for the signal events are shown. In the distribution in the top-left panel, 
we can see two bands along the x-axis at M uu — 130 GeV and 170 GeV. In the distribution in the 
top-right panel, two peaks can be found around (M TT , M uu ) = (130,170) and (170,130) [GeV]. 
In the lower panels we show the distributions of the sum of the background processes for the 
reference. The background events after solving the tau lepton momenta mostly have M TT ~ mz- 
Thus by the m^-window cut on M TT , more rejections of the background are expected. This two 
dimensional invariant mass distribution gives an evidence for the pair production of neutral Higgs 
bosons which dominantly decay into the pair of tau leptons and occasionally into the pair of 
muons directly. Accurate mass measurement is also possible, when we have sufficient statistics. 

Similar analysis to the 2/x2r^ channel can be performed for the other four-lepton channels with 
two or more muons, which are 2/zlTftle, 3/xlr^, 4/x, 3/xle and 2/x2e channels. In the channels where 
there are three or more muons, combinatorial complexity may be avoided by choosing the pair of 
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FIG. 3: Invariant mass distributions of the tau-jet pair M ThTh (left) after the selection cuts up to the 



opposite-charge muon pair in TABLE II], and those of the tau-lepton pair reconstructed by using the 
collincar approximation M TT (right) in the 2fi2rh channel. Only the events with a solution of the collinear 
approximation in < £1,2 < 1.1 ranges are used for the latter distribution. Signal events plotted in the solid 
black line are generated with m H — 130 GeV and m A = 170 GeV. Contributions from the VV background 
are also shown in the dashed red line. 



Lepton channels 


2fi2r h 
s/b (S) 


2/ilT>jle 
s/b (S) 


s/b (S) 


4 M 

s/b (S) 


3/ile 
s/b (S) 


2,u2e 
s/b (S) 


Pre-selection 
p^ h > 40 GeV 

X T > 30 GeV 
< 50 GeV 

i4 cp > 250 GeV 

(m z ) MM ± 10 GeV 


98.7/16507. (0.8) 
42.8/618. (1.7) 
34.7/ 59.7 (4.5) 
25.2/17.1 (5.1) 
22.7/13.3 (5.1) 
19.8/6.4 (5.8) 


35.7/154. (2.8) 
23.6/46.2 (3.2) 
19.2/25.3 (3.3) 

(-) 

15.7/18.4 (3.3) 
13.7/2.3 (5.9) 


14.3/162. (1.1) 
9.4/46.8 (1.3) 
7.7/27.5 (1.4) 

(-) 

6.5/20.2 (1.4) 
5.4/5.1 (2.1) 


0.9/401. (0.0) 

(-) 

0.7/4.2 (0.3) 

(-) 

0.6/2.1 (0.4) 
0.5/0.5 (0.6) 


2.6/8.9 (0.8) 
. (_) 

2.0/5.8 (0.8) 
. (_) 

1.5/2.8 (0.8) 
1.3/0.6 (1.3) 


3.3/647. (0.1) 

(-) 

2.5/5.7 (1.0) 
(-) 

1.7/3.0 (0.9) 
1.3/0.3 (1.6) 


O.S. muons 
< 21,2 < 1-1 
{m z ) TT ± 20 GeV 


13.8/6.3 (4.4) 
3.0/1.6 (1.9) 
2.8/0.8 (2.3) 


9.5/2.0 (4.6) 
3.0/1.0 (2.3) 
2.8/0.3 (3.0) 


(-) 

2.0/1.3 (1.5) 
1.8/0.7 (1.7) 


(-) 

0.4/0.2 (0.7) 
0.3/0.1 (0.8) 


- (") 
0.7/0.3 (1.0) 
0.6/0.2 (1.0) 


0.9/0.3 (1.3) 
0.4/0.2 (0.8) 
0.4/0.0 (-) 



TABLE IV: Summary of event rejections for the four-lepton channels with two or more muons. The signal- 
to-background ratio s/b and the significance S for L = 100 fb _1 are listed. The cut on H^* 1 is not imposed 
for the channels with one or less tau-jet. A pairing rule of dimuons for the channels with three or more 
muons is explained in the text. In the 2/i2e channel, an additional cut of \M ee — m z \ > 10 GeV is applied 
at the cut of {m^)^ ± 10 GeV. 



opposite-sign muons which gives the largest transverse momentum of the muon pair. The signal- 
to-background ratio and the expected significance for these channels by applying the selection cuts 
similar to those for the 2/x2r^ channel are summarized in TABLE [Fv| for L = 100 £b _1 . The cut 
variables are not optimized for each channel. The cut on is not imposed for the channels with 
one or less tau-jet, since for these channels ti process are negligible from the pre-selection level. 
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FIG. 4: Two dimensional plots for the invariant mass distributions M w vs. M ThTh (upper-left) and M w 
vs. M TT (upper-right) in the signal processes pp — > HA, where m H = 130 GeV and m A = 170 GeV are 
taken. The distributions for the background events are also given in the lower panels. The normalization of 
the box-plot in each panel is set separately. 



The collinear approximation is applied for the two leptons other than the two opposite-sign muons. 
The large s/b and S can be observed for the 2^1-rftle channel as well, where the tt background is 
less significant than for the 2fi2rh channel. 

By the requirement of opposite-sign charges in the dimuon, the background processes hardly 
lose events, while the signal process reduces one third of events which is consistent with a naive 
expectation by counting combinatory. It means, in the other word, that the requirement of same- 
sign charges in the dimuon extracts the signal events with only small background contribution 
mainly from the VV process. This is also useful to find the evidence on top of the SM process 
in the various channels with two muons or also two electrons. Those channels also appear in the 
latter. 
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FIG. 5: Distributions of kinematical variables for the 4r/,. channel. Line descriptions follow those in FIG. [I] 
2. Four-lepton channels with two or more tau-jets 

Due to the dominant branching ratio into tau leptons, a large number of events is expected 
in the channels with high tau-jet multiplicity, thus an excess beyond the SM expectation may be 
detectable in these channels. On the other hand, four-lepton channels without a primary muon 
pair from the decays of neutral Higgs bosons are hard to be kinematically reconstructed. 

First, we consider the At^ channel. Since this channel contains only hadronic objects but no 
leptonic objects in an event, the triggering efficiency may not be as good as those in the channels 
with muons. However, the requirement of the high-p T tau-jet shall stabilize the abstraction 
of the events in this channel, and we expect at least relative information would be available 
such as kinematical distributions and the signal-to-background ratio. Then, we present results 
and comments for the other four-lepton channels with two or more tau-jets. In the end of this 
subsection, we also comment on the remaining four-lepton channels. 



After the pre-selection with the requirement of the vanishing charge sum of the four hadronic 
tau-jets, the expected numbers of events are 324, 147, 797 and 5105 for HA, VV, ti and F+jets 
production processes, respectively, for L = 100 fb _1 . The dominant background contributions 
come from Z+jets production followed by the Z — > t + t~ decay and the hadronic decays of the 
tau leptons with two more tau-jets which are mis-identified from the ISR jets. 
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Arh event analysis 


HA 




VV 


tt 


V+jcts 


s/b 


S (100 fir 1 ) 


Pre-selection 


324. 


52.8 


147. 


797. 


5105. 


0.1 


4.7 


Pr h > 40 GeV 


67.2 


4.9 


2.0 


14.7 


21.7 


1.9 


9.4 


X T > 30 GeV 


48.6 


4.4 


1.1 


7.6 


10.4 


2.8 


9.3 


< 50 GeV 


34.2 


3.4 


0.5 


0.8 


8.2 


3.9 


8.7 


H^ p > 350 GeV 


27.6 


2.7 


0.4 


0.5 


3.1 


7.5 


9.3 



TABLE V: Table for the background reductions in the 47), channel. 



In FIG. ||, we show distributions of the transverse momentum pJfi of the tau-jet, that pTp Th for 
the tau-jet pair, the missing transverse energy Xt an d the scalar sum of the hadronic and leptonic 
transverse momentum, and H^ p , respectively. Distributions for all the signal and background 
processes are normalized to be unity. Line descriptions follow those in FIG. |l[ In order to reject 
backgrounds, the following selection cuts are applied; 

p% > 40 GeV, (10a) 

Xt > 30 GeV, (10b) 

H!? < 50 GeV, (10c) 

H l ° p > 350 GeV. (lOd) 

Reductions of the signal and background events are summarized at each cut in TABLE |V[ For 
L = 100 fb _1 , an expected significance S is about 5 at the pre-selection level, and reaches to 
around 9 after the selection cuts. The signal-to-background ratio s/b rises to about 8 at the end of 
all the selection cuts. Thus, finding the evidence in this channel seems very promising. It may be 
even possible at more earlier stage at the LHC. The more tight cut for the p T of the tau-jets can 
enhance the signal-to-background ratio, while slightly reduces the significance. We note that we 
have not estimated the background contributions from pure QCD processes as well as the V+jet 
production process followed by hadronic decays of weak gauge bosons, which could survive after 
the selection cuts if all the four tau-jets are mis-identified. We expect that a severe 0t cut should 
sufficiently reduce the number of events for these processes. 



Since in the Ath events all the four hadronic decays of tau leptons could yield missing momenta, 
the collinear approximation analysis used in the 2^2r^ channel cannot be applied, and the masses 
of the Higgs bosons cannot be reconstructed kinematically. However, as we will see below, distri- 
butions of the invariant mass of the tau-jet pair are useful to obtain the information of the mass of 
the Higgs bosons. To resolve combinatorial complexity from the four tau-jets, one invariant mass 
is constructed from the tau-jet pair with opposite-charges which gives the largest p^ Th , and the 
other is constructed from the remaining tau-jet pair. Then, the invariant masses of the two tau-jet 
pairs, M ThTh and m ThTh , are defined by M ThTh > m ThTh . 

In FIG. [f| distributions of M ThTh (left) and m ThTh (right) are shown for the signal events from 
HA production, while the background contributions are omitted because they are negligible. In the 
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FIG. 6: Distributions of the invariant mass of the opposite-charge tau-jet pair, M ThTh (left) and m ThTh 
(right), where M ThTh > m ThTh in the Arh channel. Signal events are generated with m H = 130 GeV and 
m A = 170 GeV. 

M ThTh distribution in the left panel, we can see an endpoint of the distribution at M TfiTh ~ 170 GeV 
which corresponds to m A = 170 GeV. Another mass information for m H = 130 GeV can be seen as 
an endpoint in the m ThTh distribution in the right panel or may be found as a bump in the M T , T , 
distribution. Thus, if sufficient numbers of events are obtained by accumulating the integrated 
luminosity such like a few thousand fb^ 1 at the LHC, it should be possible to extract the mass of 
the Higgs bosons by these distributions. In FIG. fj], we also present the two dimensional distribution 
of the invariant masses M ThTh and m ThTh for the signal events in the 4t/i channel after the selection 
cuts. 



Similar analysis to the 4rh channel can be performed for the other four-lepton channels 
with two or more hadronic tau-jets; such as 3r/jl^, 3r/jle, 27/jl/ile and Ir^le. channels. Due 
to the expected large number of events, many of them can give large S with large s/b by the 
similar selection cuts for L = 100 fb _1 . Measurements of these channels would provide a test 
of the signatures of tau leptons in the extensive channels, and thus they are useful to probe the 
production of the tau-lepton-specific Higgs bosons. The endpoint study such like in FIG. || may 
not be suitable for the invariant mass distributions including ji or e, because the momentum 
fraction z for the leptonic decays of tau leptons distributes mainly in a small z region fl5|| . The 
rest of the four-lepton channels are lr/,,l/i2e, lr^e, l/x3e and 4e. However, except the lr^l^e 
channel, the expected numbers of events are too small; a few events or less for L = 100 fb _1 . The 
signal-to-background ratio and the estimated significance by applying the selection cuts for these 
channels are summarized in TABLE VI and VII assuming the integrated luminosity of 100 fb _1 . 



Before closing the subsection for the four-lepton channels, we point out an interesting com- 
parison between the 2/i2r/j and 2rf l 2e channels. Theoretically, the difference of the two channels 
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FIG. 7: Two dimensional distributions of the two opposite-sign tau-jet pair (M ThTh vs m ThTh ) in the Ath 
channel after the selection cuts for the signal processes pp — > HA with m H = 130 GeV and m A = 170 GeV. 



arises from the direct decay of Higgs bosons into dimuons through the Yukawa couplings of Higgs 
bosons to muons. On the other hand, those to electrons are negligible, since /ml ~ 2.3 • 1(T 5 . 
The invariant mass distribution of the dielectron in the 2T/ l 2e channels is nothing but that of the 
dimuon in the 2t/j2/x channel except the absence of the two resonant peaks seen at M aa ~ m H and 
Man, — m A- The common part of the distributions between the 2/^2^ and 2rh2e channels should 
be originated from the 4r decay, while the excess in the 2/i2r^ channel relative to the 2r/j2e channel 
should be originated from the decay of the Higgs bosons directly into muons. 

Experimentally, the efficiencies of finding electrons and muons are different due to the different 
isolation condition and acceptance cuts. In our simulation study, the normalization of the events 
could be corrected by naively multiplying for the number of events in the 2r/j2e channel by 1.1, 
which is about the ratio of the finding efficiency for muons to that for electrons obtained by 
comparing the numbers of events in the 3r/jl^ and Sr^le channels. 



C. Charged Higgs boson associated production 

At the hadron collider, charged Higgs bosons can be produced in association with the neutral 
Higgs bosons via the pp — >■ W ± * — >■ (jPH^ processes f22]| . In the Type-X THDM, more than 99% 
of charged Higgs bosons decays into a tau lepton and a neutrino, and 0.35% ~ m^/m^ does into 
a muon and a neutrino for tan /3 > 2 . Therefore, the characteristic signatures for this process 
would be the three-lepton channels including many tau leptons or muons accompanied with large 
missing momentum. 

There are ten kinds of the three-lepton channels. Following the analysis of the four-lepton 
channels, we first present an analysis for the three-lepton channels with two or more muons, then 
we present analysis for those with two or more tau-jets. Brief summary of the analysis for the rest 
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Lepton channels 


477, 

s/b (S) 


3r h lfi 
s/b (S) 


377, le 
s/b (S) 


277,1/ile 
s/b (S) 


2r^2e 
s/b (S) 


Pre-selection 
pTf > 40 GeV 

X T > 30 GeV 
flff* < 50 GeV 
i4 cp > 350 GeV 
(mz)ee ± 10 GeV 


377./6050. (4.8) 
72.1/38.5 (9.5) 
53.0/19.0 (9.3) 
37.6/9.6 (8.7) 
30.3/4.0 (9.3) 
(-) 


302./4208. (4.6) 
87.2/70.2 (8.9) 
69.3/54.6 (8.0) 
49.0/17.4 (8.9) 
34.5/8.4 (8.4) 
(-) 


278./3883. (4.4) 
80.2/72.2 (8.2) 
63.4/53.8 (7.5) 
44.9/23.0 (7.6) 
31.4/10.9 (7.2) 
(-) 


166./917. (5.3) 
71.7/67.5 (7.6) 
58.0/58.6 (6.7) 
41.7/13.7 (8.5) 
24.2/3.8 (8.0) 
(-) 


74.4/13202. (0.6) 
32.4/479. (1.5) 
26.3/38.6 (3.8) 
18.7/16.0 (4.0) 
10.7/8.2 (3.2) 
9.3/2.5 (4.2) 



TABLE VI: Summary of event rejections by the selection cuts for the four-lepton channels with two or more 
tau-jets. 



Lepton channels 


lr/jl^e 
s/b (S) 


lr/j3e 
s/b (S) 


l/i3e 
s/b (S) 


4e 

s/b (S) 


Pre-selection 
> 40 GeV 
X T > 30 GeV 
(m z ) ee ± 10 GeV 


29.2/132. (2.5) 
19.3/38.6 (2.9) 
15.5/22.1 (3.0) 
13.6/2.4 (5.8) 


8.7/120. (0.8) 
5.6/34.2 (0.9) 
4.6/19.2 (1.0) 
4.0/6.5 (1.4) 


1.7/7.6 (0.6) 

- (") 
1.2/3.4 (0.6) 
1.1/1.2 (0.9) 


0.4/268. (0.0) 
(-) 

0.3/2.6 (0.2) 
0.2/0.7 (0.2) 



TABLE VII: Summary of event rejections by the selection cuts for the rest of the four-lepton channels. 
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FIG. 8: The dimuon invariant mass distribution for the 2/zlr^ channel after the pre-selection (left) and 



after the selection cuts before the cut on the m z -window in TABLE VIII (right). For the signal process, 
m H = 130 GeV, m A — 170 GeV and m H ± = 150 GeV are taken. 



of the three-lepton channels is presented after a while. 



1. Three-lepton channels with two or more muons 



In the 2/ilr/j channel, the two muons can be originated from the decay of H or A, or one muon 
comes directly from and the other comes secondarily through the leptonic decay of the tau 
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event analysis 


HA 




VV 


tt 




s/b 


S (100 ftr 1 ) 


Pre-selection 


103. 


341. 


10380. 


22683. 


2297490. 


10- 4 


0.3 


> 40 GeV 


71.4 


230. 


3147. 


6690. 


455178. 


10~ 3 


0.4 


X T > 80 GeV 


25.0 


110. 


244. 


2659. 


260. 


10~ 2 


2.4 


< 30 GeV 


10.3 


62.8 


122. 


96.5 


102. 


0.2 


3.9 


> 350 GeV 


5.9 


33.8 


50.8 


26.0 


36.0 


0.4 


3.5 


(m z ) MM ± 10 GeV 


5.4 


30.8 


17.4 


23.0 


34.0 


0.5 


3.9 


O.S. muons 


4.1 


21.3 


16.3 


23.0 


34.0 


0.3 


2.8 


> 80 GeV 


2.7 


14.4 


10.5 


12.4 


18.5 


0.4 


2.5 


(m^oWilO GeV 


0.7 


5.4 


1.9 


2.7 


1.0 


1.1 


2.3 



TABLE VIII: Table for background reductions in the 2/ilr/j channel. 



leptons. In the former case, the two muons must have opposite charges, and the distribution of 
their invariant mass shows resonance peaks at the mass of neutral Higgs bosons. 

In the left panel in FIG. [8|, we plot the invariant mass distribution of the muon pair for the 
channel after the pre-selection. The distributions are scaled to the expected number of the 
events for L = 100 fb _1 for each signal and background process. The expected number of events 
from the signal (fpH^ (and HA) process is four orders of magnitude smaller than that for the 
background processes, where the dominant contributions come from the Z+jets process. 

To extract the signal events, we consider the selection cuts on this channel. In FIG. p|, we plot 
distributions of the transverse momentum pj, for a muon, that pTp for a tau-jet and that pJ£ Th for a 
tau-jet pair, the missing transverse energy J?t, and the scalar sum and H^P of hadronic and 
leptonic transverse momentum, respectively, in the events after the pre-selection. Looking 

at the distributions in FIG. ^, we perform the following selection cuts to extract signal events in 
the 2/zItTj channels; 

p T T h > 40 GeV, (11a) 
Xt > 80 GeV, (lib) 

< 30 GeV, (11c) 
flJP > 350 GeV, (lid) 

- m z \ > 10 GeV. (lie) 



Reductions of the number of events are summarized in each step in TABLE VIII. After the m z - 
window cut on the distribution, the number of events for signal processes is ~ 36, while the 
number of background events is expected to be ~ 74. Although the s/b ratio is less than one, the 
significance S reaches to ~ 4 for L = 100 fb _1 . Therefore, if a reliable background estimation is 
achieved, the excess of the observed events could be the signature for the production of the neutral 
and charged Higgs bosons. 

In the right panel in FIG. ||, we plot the distribution after the selection cuts up to the 
cut on i?y P . After the selection cuts, the M w distribution in the y+jets background events tends 
to drop quickly before ~ 50 GeV. This is due to the requirement of the large 0t, thus the 
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FIG. 9: Distributions of kinematical variables for the 2[i\Th channel. Line descriptions follow those in 
FIG. @. 

events with a muon pair which directly comes from the decay of Z bosons are mostly rejected, 
but the events with a muon pair from the Z -4 t + t~ decay followed by the muonic decays of the 
tau leptons remain with one mis-identified tau-jet in the 2/zlr/j channel. The VV backgrounds 
mostly come from WZ production followed by Z — > ^ + /J,~ and W ± — >■ decays. Thus, the 
M MA1 distribution has a peak at M MM ~ m z , and the VV background events are reduced by the 
m^-window cut significantly. The signal events from (fpH^ production could show up sharp peaks 
at M w ~ m H and ~ m A on top of the background contributions. The resolution would 

be sufficient to determine the mass of the Higgs bosons, since it is constructed from the muon 
momenta. 

In order to measure the mass of the charged Higgs boson, we propose an analysis using the 
transverse mass variable. First, we further purify the events by requiring opposite-sign in the 
muon pair to extract the dimuon which comes from the decay of the neutral Higgs boson. For the 
same purpose, the cut on the transverse momentum of the dimuon is applied. This cut reduces 
the significance S slightly, however, it rather enhances the signal-to-background ratio s/b. Such 
purification of the signal events enables us to see the delicate structure of the transverse mass 
distribution for the signal events. The transverse mass of the tau-jet and the missing transverse 
momentum is constructed by 

M T = J2p§ -Xt • (1 - cosA(/>), (12) 
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FIG. 10: Two dimensional distributions of the dimuon invariant mass and the transverse mass of the tau-jet 
and the missing transverse momentum for the signal process (left) and the sum of the background processes 
(right) in the 2/ulr/, channel after the all selection cuts up to pt^ > 80 GeV in TABLE VIII . The signal 
process is simulated with m H = 130 GeV, m A — 170 GeV and m H ± = 150 GeV. 



where Acj) is the difference of azimuthal angles of the tau-jet and the missing transverse momentum. 

In FIG. [l(], two dimensional distributions of M^ vs. Mt are shown separately for the signal 
events (left) and the background events (right). The signal distribution shows two bands in < 
Mt < 150 GeV at M w ~ 130 GeV and 170 GeV. Clearly, these bands indicate the signature 
of (fpH^ production followed by (jP — > [i' ~fj,~ and — > t^v with the hadronic decay of the 
tau lepton. The transverse mass distribution in each band shows an edge at the mass of the 



charged Higgs boson. In FIG. 11, we show the transverse mass distributions for the mass window 
of |M MM — m^o | < 10 GeV for <fi° = H (left) and cjp = A (right). Both distributions have an 
edge at Mt — 150 GeV, which corresponds to the mass of the charged Higgs boson. However the 
expected number of events for L = 100 fb _1 is not large enough. If a sufficient number of events 
is accumulated, it should be possible to measure the mass of the charged Higgs boson from the 
position of the edge in the transverse mass distribution around M™ ~ m H and m A . 

Similar channels to the 2/xlr^ channel are 3/j and 2/^le. We summarize the selection cuts for 



these process in TABLE IX 



2. Three-lepton channels with two or more tau-jets 



Then, we study the three-lepton channels with two or more tau-jets. Similarly to the four-lepton 
channels with two or more tau-jets, these channels are expected to have large number of events due 
to the dominant branching ratios into tau leptons. On the other hand, due to the many sources of 
the missing momentum, namely neutrinos from the decay of the charged Higgs boson and the tau 
leptons, the reconstruction of the event kinematics is limited. 

Here, we present the analysis on the 3r^ channel and study how the mass information can be 
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FIG. 11: Distributions of the transverse mass of the tau-jet and the missing transverse momentum with 
120 GeV < M MM < 140 GeV (left) and 160 GeV < M w < 180 GeV (right) for the 2//lr h channel after the 
selection cuts. The signal processes is simulated with m H — 130 GeV, m A — 170 GeV and m H± = 150 
GeV. 



Lepton channels 


2/xIt^ 
s/b (S) 


3/i 

s/b (S) 


2/tle 
s/b (S) 


Pre-selection 
Pr h > 40 GeV 
X T > 80 GeV 
< 30 GeV 
H^ p > 350 GeV 
(ro z ) M/1 ± 10 GeV 


445./2330550. (0.3) 
301./465015. (0.4) 
135./3163. (2.4) 
73.1/320. (3.9) 
39.7/113. (3.5) 
36.2/74.4 (3.9) 


34.1/4214. (0.5) 
(-) 

13.9/416. (0.7) 

(-) 

7.5/122. (0.7) 
7.0/67.8 (0.8) 


81.8/3880. (1.3) 
(-) 

30.6/385. (1.5) 
(-) 

14.3/115. (1.3) 
13.2/11.9 (3.3) 


O.S. muons 
> 80 GeV 
(m^ilO GeV 


25.3/73.3 (2.8) 
17.1/41.5 (2.5) 
6.1/5.4 (2.3) 


(-) 
(-) 

3.0/7.7 (1.0) 


9.2/10.0 (2.6) 

(-) 

4.0/1.0 (2.9) 



TABLE IX: Summary of event rejection by the selection cuts for the three-lepton channels with two or more 
muons. 



extracted in this channel. In FIG. 12, we show distributions of kinematical variables, , p^ Th , 
$Ti an d H^ p . The numbers of the signal and background events are normalized to be unity. 
Taking into account these distributions, we apply the following selection cuts for the signal and 
background events; 

p% > 60 GeV, (13a) 
Xr > 80 GeV, (13b) 
< 50 GeV, (13c) 
> 350 GeV. (13d) 
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FIG. 12: Distributions of kinematical variables for the ?>Th channel. Line descriptions follow those in FIG. 



3r/j event analysis 


HA 


o i/± 


VV 


tt 


y+jets 


s/b 


S (100 ft- 1 ) 


Pre-selection 


710.5 


2476. 


7069. 


42053. 


777043. 


10~ 3 


3.5 


> 60 GeV 


88.4 


299. 


34.2 


341. 


1237. 


0.2 


9.3 


X T > 80 GeV 


29.3 


165. 


9.9 


89.9 


272. 


0.5 


9.3 


ffjf < 50 GeV 


15.8 


114. 


4.0 


4.8 


123. 


1.0 


9.9 


H l ° p > 350 GeV 


15.0 


106. 


3.4 


3.8 


108. 


1.0 


9.9 



TABLE X: Table for background reductions in the 2>Th channel. 



The summary of the results of event reductions is given in TABLE [X|. After the all cuts, more than 
a hundred event is expected for the signal processes, while the background events are estimated 
to be around 120 which is dominated by the V+jets process for L = 100 fb -1 . The expected 
reach of the significance is close to 10 for L = 100 fb . It may indicate the possibility of an 
earlier evidence of the signal in this channel. The larger H^ p cut can be applied to enhance the 
signal-to-background ratio, which however tends to reduce the significance. 

In the left panel in FIG. 13, we plot the invariant mass distributions of the tau-jet pair with 
opposite charges after the selection cuts for the signal and background processes, where the pair 
of the two tau-jets is chosen such that it has opposite charges and gives the largest transverse 
momentum p^ Th . The M ThTh distribution shows an endpoint behavior at M ThTh ~ 170 GeV which 
corresponds to the mass of the heavier neutral Higgs boson. Another sharp edge may be able to 
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FIG. 13: Distributions of the invariant mass of the tau-jet pair (left) and the transverse mass of the remaining 
tau-jet and the missing transverse momentum in the 3t/j channel after the selection cuts. The pairing rule 
of the tau-jets is explained in the text. 



be pinned down at M ThTh ~ 130 GeV the mass of the lighter neutral Higgs boson. In the right 
panel, distributions of the transverse mass Mt constructed from the remaining tau-jet and the 
missing transverse momentum using Eq. ( |l2| ) are shown for the signal and background events after 
the selection cuts in the 3rt channel. The Mt distribution for the signal events vanishes around 
Mt — 150 GeV, however the slope is, in contrast to the 2//1t^ channel, rather gentle but not seen 
like an endpoint. This is because there are many sources of the missing momentum in this channel. 

In the left panel in FIG. 14, we plot the two dimensional distribution of the invariant mass of 
the tau-jet pair and the transverse mass in the 3^ channel after the selection cuts. In the right 
panel, the dominant F+jets background events are plotted. The V+jets background events widely 
distribute in the two dimensional space of M ThTh and Mt- On the other hand, the signal events 
distribute compactly in the region of 50 GeV < M ThTfi < max (m#, m^) and < Mt < rn H ±- 
Therefore, if the sufficient number of events is accumulated, we may have an indication of the 
masses of the charged Higgs boson and the neutral Higgs bosons by catching the border of the 
signal-event-like distribution. 



Similar channels to the 3th channel are 2r/il// and 2r/jle. We summarize the selection cuts 
for these processes in TABLE XL The three- lepton channels with one or less muon or electron 
suffer large background contributions from the W+jets process, thus the s/b ratio is not so large. 
However due to the large number of the signal events, the significance S can be large enough. 
There are also lr^l/xle, lr^,2e, l/x2e and 3e channels. For the reference, we list the summary of 
reductions of events in TABLE XII. 
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FIG. 14: Two dimensional distributions of the invariant mass of the tau-jet pair versus the transverse mass 
in the ?>Th channel after the selection cuts. The left panel is for the signal process, and the right panel is for 
the background processes. The normalization of the two panels is adjusted. 



Lepton channels 


3r/j 

s/b (S) 


2r h ln 
s/b (S) 


27), le 
s/b (S) 


Pre-selection 
Vt > 60 GeV 
X T > 80 GeV 
< 50 GeV 
H l ° p > 350 GeV 


3186./826165. (3.5) 
387./1612. (9.3) 
194./372. (9.3) 
129./132. (9.9) 
121./116. (9.9) 


2001./622055. (2.5) 
494./7346. (5.7) 
244./1940. (5.4) 
158./761. (5.5) 
124./434. (5.7) 


1794./571334. (2.4) 
443./7368. (5.1) 
216./1805. (5.0) 
142./738. (5.1) 
112./437. (5.2) 



TABLE XI: Summary of event rejection by the selection cuts for the three-lepton channels with two or more 
tau-jets. 



IV. SUMMARY AND DISCUSSIONS 



In this section, we present a summary and discussions of our simulation studies, and give 
discussions for the neutral Higgs pair production at electron-positron linear colliders. Finally, we 
conclude the paper. 

We have studied multi-r signatures at the LHC in the Type-X THDM. Assuming tan f3 > 3, 
more than 99% of extra Higgs bosons H, A and decays into tau leptons. These Higgs bosons 
can be pair produced in qq — > HA and qq' — > (j)°H processes where (fP = H, A, and they 
produce characteristic multi-r final states at the LHC. We have focused on the 2/x2r and 4r states 
for the pair production of the neutral Higgs bosons, and on the 2pAT\v and 3tIv states for the 
charged Higgs boson associated production followed by the decays of tau leptons into leptons or 
hadrons. These result in the various four- and three-lepton channels which include tau-jets. We 
have performed the Monte-Carlo simulation for the signal processes as well as the major background 
processes, namely VV, ti and V+jets processes at the LHC where V = W, Z. The realistic tau-jet 
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Lepton channels 


lT^l/^le 
5/6 (S) 


lr h 2e 
s/b (S) 


l/i2e 
s/b (S) 


3e 

s/b (S) 


Pre-selection 
Vt > 40 GeV 
X T > 80 GeV 
< 30 GeV 
> 350 GeV 
(m z ) ee ± 10 GeV 


738./65150. (2.9) 
498./17718. (3.7) 
302./4971. (4.2) 
113./455. (5.1) 
58.5/128. (4.8) 
(-) 


333./1883890. (0.2) 
226./373858. (0.4) 
97.3/2456. (2.0) 
52.8/259. (3.2) 
27.7/94.4 (2.7) 
24.2/62.8 (2.9) 


68.8/3431. (1.2) 
(-) 

22.7/343. (1.2) 
(-) 

9.1/99.5 (0.9) 
7.9/11.2 (2.2) 


20.0/3107. (0.4) 
(-) 

6.7/302. (0.4) 
(-) 

2.7/90.7 (0.3) 
2.2/39.9 (0.3) 



TABLE XII: Summary of event rejection by the selection cuts for the rest of three-lepton channels. 



tagging procedure is simulated for the final state hadrons. For leptons, the standard acceptance 
cuts and isolation conditions are considered. Then we have applied the selection cuts to reduce 
the background contributions for each channel. We here summarize the analysis for the four- and 
three-lepton channels to probe the pair production of the neutral Higgs bosons and the charged 
Higgs boson associated production, respectively, in order. 

• For the four- lepton channels with two or more muons, the mass of the neutral Higgs bosons 
can be measured as the peaks in the M w distribution. For these events where the dimuon 
comes from the direct decay of the Higgs boson, the tau lepton momenta can be recon- 
structable by using the collinear approximation for the rest two leptons out of the four. 
The two dimensional distribution of the Muu an d the reconstructed M rr , such like shown 
in FIG. [|, would be a direct test of the neutral Higgs bosons pair production in the Type-X 
THDM. However, the numbers of the signal events are relatively small due to the small 
branching ratio of the neutral Higgs bosons into dimuons. The signal-to-background ratio 
is rather large, thus a clean extraction is possible, if the number of events is sufficient. The 
information of the branching ratio of the neutral Higgs boson into dimuons can also be 
extracted. 

The detailed analysis is performed for the 2fj,2rh channel, and the summary of the selection 



cuts for the similar channels is found in TABLEs [II and IV 



• For the four-lepton channels with two or more tau-jets, a large number of signal events are 
expected due to the dominant branching ratio of the neutral Higgs bosons into tau leptons. 
After the selection cuts, both the signal-to-background ratios s/b and the significances S 
are large, so that the excess can be relatively easily found. On the other hand, the mass 
determination would only be possible by looking at the kinematical distributions. Searches 
for the endpoints or bumps in the M ThTh distributions, shown in FIGs. [6| and [7| can indicate the 
mass of the neutral Higgs bosons, provided the sufficient number of events are accumulated. 

The analysis for the 4r^ channel is presented in detail, and the selection cuts for the similar 
channels are summarized in TABLEs |V] and |VT| . For completeness, a summary of the selection 



cuts for the rest of the four-lepton channels is presented in TABLE VII 
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For the three-lepton channels with two or more muons, the mass of the neutral Higgs bosons 
can be measured by looking for the peaks in the Man distribution. The mass of the charged 
Higgs boson can be measured as the edge of the distributions of the transverse mass Mt 
which is constructed from the rest of the three-leptons and the missing transverse momentum. 
By focusing on the events with Mnu ~ m H or m A , an edge can be enhanced such as shown in 



FIG. 11, which corresponds to the charged Higgs boson mass. Due to the large background 
contribution from the y+jets process, the s/b ratio is small. 

The analysis for the 2/ilr^ channel is presented in detail, and the selection cuts for the the 
similar channels are summarized in TABLEs VIH] and IX. 



For the three-lepton channels with two or more tau-jets, a large number of signal events is 
expected due to the dominant branching ratio of Higgs bosons into tau leptons. However 
due to the large background contribution from the V^+jets process, the s/b ratio is rather 
small even after the selection cuts. Using the two dimensional distributions of M ThTh and the 
transverse mass Mt, as shown in FIG. [l4|, the information of the mass of the Higgs bosons 
can be obtained. 

The analysis for the 3t/i channel is presented in detail, and the selection cuts for the similar 
channels are summarized in TABLEs |X| and pOl. For completeness, a summary of the selection 



cuts for the rest of the three-lepton channels is presented in TABLE KII 



We have shown that the characteristic excess from the production and the decay of the 
lepton-specific Higgs bosons can be observed in the various four- and three-lepton channels on top 
of the SM contributions at the LHC. However, on the other hand, the determination of the mass of 
the Higgs bosons seems not to be feasible with the integrated luminosity of L = 100 fb , although 
we have shown the various methods in various channels. The measurement of the branching ratio 
of the neutral Higgs bosons into dimuons which we have examined in the 2[i2Ty l channel, also 
requires the number of events more than the expected for L = 100 fb -1 . These measurement 



should be achieved with the running of the luminosity upgraded LHC in future [46]. With the 
luminosity of an order of thousand inverse fb, it is promising to measure the mass of the neutral 
and charged Higgs bosons as well as the branching ratio of the neutral Higgs boson into dimuons. 

We here briefly comment on the search potential for the lepton-specific Higgs bosons at the 
International Linear Collider (ILC) or the Compact Linear Collider (CLIC). The neutral Higgs 
bosons can be pair produced via the e + e~ — > HA process, and their decay produces four-r 
states predominantly. Through the leptonic and hadronic decay of the tau leptons, the signatures 
are four-lepton channels including tau-jets, as we have studied for the hadron colliders. At 
e + e~ colliders, four momenta of the four tau leptons can be solved by applying the collinear 
approximation to all the four decay products of the tau leptons |25|, H7j, because the missing four 



momentum can be reconstructed by the energy momentum conservation. In Ref. [48|, it is shown 
that clear peaks can be observed in the one dimensional and two dimensional distributions for the 
invariant masses of the reconstructed tau lepton pairs. 
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In conclusion, we have studied the multi-r signatures at the LHC in the THDM with lepton- 
specific Yukawa interactions. At the LHC, such Higgs bosons can be pair produced via the weak 
boson mediated processes, and predominantly they decay into multi-r states. Although the event 
reconstruction is limited due to the missing energies from the decays of the tau leptons, it is possible 
to discriminate the signal events from the background events by requiring the high multiplicity of 
the leptons and the tau-jets with appropriate kinematical cuts in various four and three-lepton 
channels. The mass of the Higgs bosons can be measured by using the invariant mass distributions 
of the dimuon which comes from the direct decay of the neutral Higgs bosons as well as the 
endpoint analysis for the distributions of the invariant mass of the tau-jets and the transverse 
mass, if sufficiently large integrated luminosity is accumulated. We have found that measurements 
of multi-r final states at hadron colliders can be a useful probe of the properties of the Higgs bosons 
in the lepton-specific models. 
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Appendix A: Tau-jet tagging efficiency 

In this appendix, we present a rough estimation of the tau-jet tagging efficiency and the mis- 
tagging probabilities for the non-tau jet in our tau-jet tagging procedure defined in the main text. 
We neglect the r/ and p T dependences of the efficiency as well as the mis-tagging probability. This 
information may be useful to estimate the variation of our results by changing the conditions of 
tau-jet tagging, and also to compare with the actual experimental situation at the LHC. 

The tau-jet tagging efficiency is estimated from the signal events in 4r^ and 3r/j channels from 
HA production and (fpH^ production, respectively. In the 4r^ channel, the number of events 
after the pre-selection for the signal events is ~ 320 for the integrated luminosity of 100 fb _1 (see 
TABLE |V|). On the other hand, the expected number of the 4r final-state is 50 [fb] x 100 [fb _1 ] ~ 
5000. From this, the probability of finding the tau-jet per decay of one tau lepton is estimated 
to be (320/5000) 1 / 4 ~ 0.50. This can be interpreted as the product of the branching ratio of the 
hadronic decay of tau leptons, the acceptance ratio for \rj\ < 2.5 and p T > 10 GeV and also the tau- 
jet tagging efficiency. By using the branching ratio of hadronic decays of tau leptons which is about 
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65% 0], the multiplication of the the acceptance ratio and tau-jet tagging efficiency is estimated 
to be ~ 0.77. Further separation of the two factors is too confusing to us. For confirmation we 
have checked that the obtained value gives roughly correct number of the signal events from (fpH^ 
production in the 3t/i channel by 200 [fb] x 100 [fb -1 ] x (0.50) 3 ~ 2500 (compare with the number 
of signal events from (fpH^ 1 production after the pre- selection in TABLE |X|). 

The mis-tagging probability for non-tau jets is estimated for the y+jets events in the 2/^2-7"^ 
and 2/ilr^ channels. The numbers of events from the y+jets production after the pre-selection 
are 2.9T0 4 and 2.3- 10 6 for the 2//2r/i and 2^1^ channels, respectively, where here the requirement 
of the zero charge-sum is not imposed in the former channel. For our y+jets events, the fraction 
of the events which exclusively contain one (two) primal jet(s) with < 2.5 and p T > 10 GeV is 
found to be 0.09 (0.01). The expected number of events for the fJ- + ^~jj final-state from the y+jets 
process is 30 [nb] x 100 [fb -1 ] x 0.33 x 0.01 ~ 10 7 , where a factor of 0.33 comes from the fraction 
of the dimuon final-state in the leptonic Z-boson decays. Thus the mis-tagging probability per jet 
is estimated to be ^2.9 • 10 4 /10 7 ~ 0.05. The same probability is estimated for the 2/ilTh channel 
as 2.3 • 10 6 /(10 9 x 0.09) ~ 0.03. Thus the mis-tagging probability per jet for the F+jets process 
should be a few percent in our simulation. 

Similarly the mis-tagging probability for the ti events can be estimated by using the dimuon 
decay of the ti events. For the ti events, one or two jets are expected to be originated from the b 
quark, and the mis-tagging probability for such jets may be different from that for the light-flavor 
jets. However, we ignore the difference of the mis-tagging probabilities by the flavor of jets. The 
number of events for the dimuon decay mode of our ti events is given as 493 [pb] x 100 [fb -1 ] x 
(0.107) 2 ~ 5.6T0 5 . For the dimuon decay channel in our ti events, an average number of the primal 
jet with \rj\ < 2.5 and p T > 10 GeV is Yl n =i r n n ~ 3.6 where r n is the fraction of events with n 
primal-jets which satisfies Yln=o r n = 1> an d a l so an average number of the combination of the two 
primal jets is Yln=2 r « (2) ~ 5.8. By removing the requirement of zero charge-sum of the four-lepton 
at the pre-selection in the 2/i2r^ channel, the number of events after the pre-selection is ~ 800, thus 
the mis-tagging probability is estimated to be -^7800/(5. 6 • 10 5 x 5.8) ~ 0.02. On the other hand, 
for the 2^1t/i channel, the mis-tagging probability is estimated to be 2.3T0 /(5.6T0 5 x 3.6) ~ 0.01. 
Thus the mis-tagging probability per jet for the ti process should be also a few percent. The smaller 
mis-tagging probability for the ti process by our rough estimation may be due to the large QCD 
activity in the ti events which could prevent the isolation of the charged tracks in the tau-jets. 

One may be able to reduce the mis-tagging probability with minor loss of tagging efficiency, 
e.g., by using only the 1-prong (or 3-prong) decay mode of tau leptons and/or by applying more 
tight tagging conditions. 
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